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SYNTHESIS OF BLADE FLUTTER VIBRATORY PATTERNS 
USING STATIONARY TRANSDUCERS 
by A. Kurkov and J. Dicus 
National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 
ABSTItACT 

Flutter frequency was determined and rotor vibratory amplitude 
and phase distributions during flutter were reconstructed from station- 
ary aerodynamic tyix? measurements. A previously reported optical 
metliod for measuring blade-tip displacement during flutter was ex- 
tended by means of digital analysis. Displacement amplitudes and 
phase angles were determined based on this method. For selected 
blades, sixictral results were also obtained from strain gage measure- 
ments. The results from these three types of measurement were com- 
pared and criiit ally evaluated. 

NOMENCLATURE 

A amplitude in rotating frame of reference 

Aj i^*' amplitude in stationary frame 

E engine order 

F Fourier transform 

f frequency 

j imaginary unit 

N number of points 


STAR Category 
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A U 

Nj i nodal diameter 

Nj total number of nodal diameters 

n numbc'r f)f blades 

t time 

ni angular separation between two sensors 

€ defined in Apix?ndix 

0 phase an^i:le in rotating frame 

A0 phase lag due to rotation 

phase corresix>nding to amplitude 
tp angular coordinate in rotating frame 

Av’ angular separation between two blades 

u) circular frequency 

rotational frequency 

INTRODUCTION 

Most commonly, turbomachinery flutter instrumentation consists 
of blade mounted strain gages. Recently, however, two alternative 
schemes | 1|* that rely on optical measurement techniques were success- 
fully implemented in a full scale engine research facility. One scheme 
produced the display of the stroboscopic image (SI) of the tip of the blade 
during flutter. The other scheme, termed Photoelectric Scanning method 
(PES), resulted in the screen display of the extremes of the blade-tip 
displacements lor each blade. This second technique was found to be 
particularly useful for flutter monitoring. Quantitative measurements, 
however, were limited to the determination of vibratory displacement 
amplitudes based on photographs of PES screen display. 

iNumbers in brackets designated References at the end of the |)a|x?r. 
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In this pa|x?r the use uf stationary sensors such as in reference [ 1| 
is explored further. Particular attention is given to quantitative analy- 
ses which until now were limited in extent. This has been achieved 
mainly through Die use of digital analysis techniques. It is anticipated 
that with further development of analysis and instrumentation techniques, 
the reliance on strain gages and the associated signal transmission 
hardware could lx* greatly reduced. 

It is shown in the first jiart of the pai>er tliat flutter can be detected 
from the data recorded by the Air-stream (i.e. aerodynamic) instru- 
mentation. A method is presented which allows construction of the 
rotor vibratory jiattern from these measurements. It relies on digital 
spectral analysis procedures incorporating Fast Fourier Fransform 
(FFT) algorithm. A method is also presented whicli allows precise 
determination of flutter frequency based on measurements provided by 
two aerodynamic transducers. 

The second part of the jiaper deals with PES data. In the present 
case, however, the displacement are obtained by means of digital analy- 
sis of the blade sensor and reference sensor signals. These signals 
were recorded on a magnetic taix? during flutter testing. Once the rela- 
tive displacements are obtained for a number of rotor revolutions, 
sijectral analysis techniques are applied to obtain vibratory amplitude 
and phase angU s. In this case the sijectral analysis is also based on 
the FFT procedure. 

Optical and aerodynamic data were available for the same test 
point so that detailed comparison of the spectral results from both 
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melhods was possible. The results exhibited similar specific features 
associated with the fact that these sensors are stationary and not blade 
mounted. Since a number of blades were also strain gaged, both sta- 
tionary and blade mounted vibratory n easurements were available for 
these blades enabling comparison of re.^ults from these two ty|)es of 
measurements. 

It should be noted that another optical metJiod for measuring of 
blade displacements during flutter has lx?en recently implemented [2]. 
Some additional information related to the use of stationary sensors 
is available in reference | 3|. Both, flutter and resonant blade vibra- 
tions are included. 

The data presented in the i)aper were obtained during the first fan- 
rotor flutter test of a preproduction turbofan engine in an altitude test 
facility. Flutter was classified as "stall flutter." It occurred in the 
neighborhood >1 70 |X?rceni of the design six?ed at high incidence and 
under elevated inlet pressure and temjx'rature conditions. Stress 

p 

levels during Mutter were about 3.5 kN cm (5 ksi). The vibratory 
mode was identified to be predominantly torsional. 

AERODYNAMIC FLUTTER-INSTRUMENTATION 

The first fan-rotor case was instrumented with a number of high- 
resixmsc statu pressure transducers. The transducers were very 
nearly Mush-mwufi'ed. Six transducers were installed in the fan case 
approximately in the direction of the rotor blade chord at the tip sjjan. 
The first transducer was aligned with the blade leading edge and the 
sixth with the blade trailing edge. Several transducers were also in- 
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stalled in the ian case curres|x>ndin(; to rotor blade midchord p<jsition 
and s|)a('ed circumferentially one blade pilch apart. 

For a limited number of lest ptjinls, the first- rotor dischari^e 
(l.e. rotor wake) was instrumented with a high>resix)nse probe in the 
form of a miniature cantilever beam (4|. Its response was proportional 
to the lift force acting on the cantilevered beam. The measurements 
with this probe were obtained at about 90 percent of the blade s|jan. 

The responses of the cantilevered beam and the static pressure trans- 
ducers were sufficiently high to capture the blade passing frequency, 
as well as its significant harmonics. The data from these measure- 
ments were recorded on a multichannel FM magnetic taix' recorder at 
304.8 cm/s (120 in s). 

SPECTRAL ANALYSIS OF AERODYNAMIC DATA 

The S|X?ctrai analysis of aerodynamic data was i^erformed on a 
digital spec tral analyzer which could be programmed to perform special 
functions in addition to the standard functions. A summary of the param 
eters used in ihe s|X‘ctral analysis and the description of the special 
procedures incorporated in the analysis arc given in the Appendix. 

Figure 1 illusi rates the p«.*wer s|)ectrum of the wake-probe trans- 
ducer signal during flutter, and Fig. 2 in the absence of flutter. The 
engine corrected cqx'raung conditions corresponding to these figures 
were the same, however, Fig. I corresponded to a higher engine inlet 
pressure. Frequencies in these figures were normalized using shaft 
rotational frequency: that is. they were expressed in engine orders. 

To facilitate interpretation of the results, Uie frequency corresponding 
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to each significant })eak was written just above that peak. Normalized 
flutter frequency based on strain gage measurements is also included 
In Fig. 1. 

The amplitudes in Figs. 1 and 2 are normalized with respect to 
the highest ix>aks. This is found to be convenient for the |)ur|X>se of 
com liar isun of S|X>ctral results based on measurements obtained with 
different ty|x?s of transducers. 

Insiiection of these figures reveals tliat discrete nonintegral order 
frequencies are characteristic of ft»*.tter and that tliey differ from flutter 
frequency by an intregral number of engine orders, within the available 
resolution accuracy. Each nonintegral order frequency corresjxjnds 
to a different travelling wave. This follows from the fact that a Np 
nodal diameter wave motion described in a frame of reference fixed 
with respect to the rotor by 

Aj cos(u;t 1= + 0^) 

transforms into 

A| COS| (a> ± u)j,Nj)t f ii>j] 

in a stationary reference frame. The upiier and lower signs in these 
expressions correspond to forward and backward moving waves, re- 
spectively. 

Flutter (>»uld also be detected from the static pressure transducer 
signals. For example. Fig. 3 presents a ixjwer spectral plot that 
corresponds to a midrhord static pressure transducer. Comparing 
Figs. 1 and 3 it is seen that there exists a reasonable agreement among 
the relative magnitudes of flutter peaks. Spectral plots corresponding 
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to other midchoitl transducers are qualitatively quite similar. At the 
leading; edge, however, sharp attenuation of high frequency ncxlal diame- 
ter waves is observed, Fig. 4. Although the trailing edge static pres- 
sure transducer was ino|X‘rative for this test jxiint, evidence from other 
test points during which this transducer was o|K?rational indicates a 
similar attenuation of the high-frequency nodal waves. Since these 
nodal waves were well reprixluced in Fig. 1, which is based on mea- 
surements in the rotor wakes at 90 jjcrcent of the blade span, this 
suggests that the attenuation phenomenon is due to the specific features 
of the flow in the immediate vicinity of the casing wall in the blade 
trailing-edge and possibly leading-edge regions. This observation has 
obvious practical significance in cases where it is desired to use static 
pressure transducers for investigation of flutter. 

AMPLITUDE AND PHASE ANGLES 

In this section the rotor vibratory pattern including both, amplitude 
and phase angles is constructed from a pressure transducer signal. A 
different test pront is selected for this purjxjse, mainly because detailed 
optical sensor data are also available for this point. For the present 
purjxjse. It IS essential to be able to determine rotor orientation a» 
time t - 0, corresponding to the start of the data sampling. The start 
of sampling was therefore, triggered using once-ijer-revolution signal 
pulse. 

The linear s|X?(tral plot (i.e.. the magnitude of the Fourier trans- 
form) corresponding to a midchord static pressure transducer is pre- 
sented in Fig. 5. Whi’.e in Figs. 1 and 3 one could identify a pre- 
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dominantly Ihrc'^ nodal diameter forward travelling wave pattern, in 
Fig. 5 no single wave pattern la dominant. In addition to the nodal* 
wave amplitudes ’’Aj) represented in Fig. 5, one also needs the corre- 
six)nding phase angles (0|). Their determination is discussed in the 
Appc'ndix. Dasically, it involves an interpolation process which be- 
comes necessary because of an inherent feature of results obtained 
with the FFT prcKedure 1 5 |. Once A.’s and 0^*s are determined, 
transformation into a rotating frame of reference and suiX‘rposition 
is accomplished using (he foilowing equation, 

y- ^ (1) 

where A and 0 represent amplitude and phase in tlie rotating frame 
of reference. Tnd </> - 0 denotes the rotor radius |X)inting to the trans- 
ducer |X)bition at time t 0. (The radius 0 = 0, approximately corre- 
sjxmded to blade number 19.1 The sign of in the above equation is 
either positive or negative deiiending on whether the wave is forward 
or backward. Tiius, referring to Fig. 5, at frequency f - 11.45, 

= 3; and at f 5.45, - *3. Ail amplitudes corresponding to 

nonintegrai frequencies marked in Fig. 5, with the exception of the 
amplitude at f - 10.447, were included in the summation. Consequent- 
ly. was 9 in equal 1 m (D. 

The results of su(jer|»sition are presented in Figs. 6 and 7. The 
amplitudes in these figures were normalized with respect to blade num- 
ber one, and the phase angles were evaluated relative to blade number 



9 


one. Bla*-e |>i8iiions were identified with an accuracy of atxjut one hali 
of the blade pitch. Large amplitude and phase angle variations evident 
In Figs. 6 and 7 are typical of the **mistunod rotor” response IGJ. 

These results will be compared later with the strain gage and optical 
displacement measurement s . 

DETERMINATION OF FREQUENCY 
The results presented in the previous section could not have been 
derived without the kn<»wledge of flutter frequency, because there would 
be no basis for selection of N^ corres|X>nding to each nonintegral en- 
gine order |X‘ak Flutter frequency may be obtained from the strain 
gage data, or it may be known within the accuracy of one engine order. 
However, frequency can also be determined precisely and without any 
previous knowledge if two aerodynamic sensors located at different 
circumferential positions around the rotor are available. It is neces- 
sary to determine the phase difference between the outputs of two trans- 
ducers at only one nonintegral engine order frequency. This phase 
difference is then equated to noN., where a is the angular separation 
between two transducers and N^ the desired nodal diameter number. 
Once one N^ is determined, the others can be determined by observing 
relative frequencies. Table 1 illustrates these phase relationships for 
two midchord transducers with ct 47.4^. Rather than calculating N^ 
in the last column fi/r iliustiaiive purposes it is more convenien*^ to 
use tlie known N^*s to calculate the phase angles. The agreement be- 
tAveen these (columti 3* arid the phase angles obtained from direct mea- 
surement (column 2) is particularly satisfactory for backward travelling 


waves . 
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OPTICAL-SENSOR SIGNAL ANALYSIS 

A brief description of the physical arrangement of different sen- 
sors employed in the PES method is included here in order to facilitate 
the explanation of the analysis procedure. The information related to 
the optical probe design is available in reference ( 1 1. In this method, 
measurements are provided by four sensors. Two optical ty]x? sensors 
were used to detect blade tip |>osltions, while two magnetic type sensors 
provided reference signals related to the shaft rotational frequency and 
the l-lade passing frequency. One of the optical sensors was positioned 
op|X>site the blade trailing edge ixisition, and tJie other opixisite the 
blade midchord. The choice for these |X)sitions, to some extent, was 
determined by the availability of an access port in the fan case. As the 
blade tips passed by each sensor, a series of pulses were generated 
and recorded on the magnetic ta|>e. Also recorded on the magnetic 
tape were the pulse signals corres|X)nding to the rotational frequency 
and the blade passing frequency. The latter were generated by a shaft 
mounted gear having the same number of teeth as the number of blades. 

The four channels of data were recorded on a multi-channel direct - 
record type tape recorder at 304.8 cm/s (120 in s). Subsequently, these 
data were recorded or «tn FM tape recorder at 304.8 cm s (120 in/s) 
playing back tJie original tajx? at 7.62 cm/s (30 in/s). The latter FM 
tape was then digitized at 9.525 cm s (3. 75 in/s) at the rate of 10 kHz 
ix?r channel. In terms of real time, the resulting digitizing rate was 
1.28 MHz. The remaining data processing was |x?rformed on a digital 
comixiter. It involved computation of relative time of occurrence of 
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different pulses. For example, in order to compute trailing edge dis> 
placements, it was necessary to calculate the time dUference between 
the occurrences of trailing edge pulse and the closest reference blade* 
passing-frequency (Xilse, for each blade and every revolution. The 
deviations about the average of these time differences are proixirtional 
to displacements during flutter. S|x>cial care had to be taken to avoid 
occasional spurious signals which were probably due to air particle* 
contamination. The particles could either emit or reflect the light. 

The results of this procedure are Illustrated in Fig. 0 where each 
unit corres{)onds to a trailing edge displacement >f about 0.27 mm 
(0.0106 in). The data frame for this figure corres|x>nds to 27 rotor 
revolutions. (Note that some of the ixiints overlap in Fig. 8.) It 
should be noted that analogous to trailing-edgc displacements, the 
midchord, and the trailing edge relative to midchord displacements 
can be evaluated. These displacements, however, are not as large, 
so that in order to achieve a reasonable accuracy, more |X)ints |x?r 
blade passage and, therefore, a higher digitizing rate than the one 
used for Fig. 8 is necessary. 

SPECTRAL ANALYSIS OF BLADE DISPLACEMENTS 

Six?ctral analysis of the displacements represented in Fig. 8 can 
be performed in two ways. One is to include displacement from all 
blades in the order they were sampled (sampling rale 38E), and the 
other is to analyze displacemenls corres|Xjnding to each blade separate- 
ly (sampling rate IE). The overall s|)ectrum, corresponding to the 
sampling rate of 38E, is given in Fig. 9. The units in this figure are 
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the same as In Fig. 8. Because of the nature of Uk. sampting process, 
only flutter, or nonintegral engine order vibrations are delected. There- 
fore, for the purjX'se of compari.«»on of results included in Figs 5 and 9, 
only nonintegral order amplitudes in Fig. 5 are relevant. These ampli- 
tudes are seen to be in approximately the same relative proporiionh in 
the two figures. 

An example of a s|)ectrum plot for a single blade corresixmding to 
the sampling rate of IE is im luded in Fig. 10. As a direct consequence 
of the sampling theorem |5|. the frequencies in this figure are aliased 
between zero and one half engine order. Frequency information in this 
figure is, therefore, more indeterminate than in Fig, 9, that is, it is 
only known by how much flutter frequency differs from a multiple of 
1, 2E, and not IE frequency, as in Fig. 9. Prec ise determination of 
flutter frequency, a.s already discussed during the analysis of static* 
pressure-transducer data, requires two sensors. Only one trailing- 
edge sensor was available in the present investigation, so that this 
procedure c juld not be re|)eated tor the optic al data; however, liased 
on sixH tral plots one can ex|xt t similar consistency of result-, as in 
Table 1. In the absence of a syctem mode of vibrat’ons, that is, when 
each blade vibrates at a frequency different from any other blade, the 
numlx?r of sensor-c should exceed two times the flutter frequency ex- 
pressed in engine orders. For the system mode, there is an up|)er 
limit on the highest nodal diameter which can be detected wiih one 

n 2 -1 


sensor: 
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whe»*e f IS the flutter frequency In engine orders. 

Referring to Fig. 9, an additional interpretation of spectral ijeaks 
In this figure is |x>ssible based on the orthogonality property of nodal 
waves. Squaring and integrating real and Imaglnery parts In equation (1) 
Ijetween the limlt^ of *11 to II, one obtains 



which relates the amplitudes in Fig. 9 to the rotor-mean squared dis- 
placement amplitudes. 

By jX'rforming spectral analysis se()arately for each blade, dis- 
placement amplitude and phase distributions around the rotor can be 
constructed. It should be noted, however, tliat in relating the phase 
angles from different blades, an account has to be taken of the fact 
that blades were not sampled simultaneously. Between any two blades 
there exists a phase lag due to the rotation. It is given by 

A0 • f 

where is the angular separation between two blades, a id I is 
flutter frequency expiessed m engine orders. The amplitude and phase 
angle distntiution-. are presented in the next section where they are also 
compared with corresixmding results from pressure and strain-gage 
measu rements. 

DISCUSSION OF RESULTS 

The normalized amplitudes corresponding to pressure, displace- 
ment, and stress measurements aie given in Fig. 11, and the relative 
phase angle.s corres|X)nding to these measurements in Fig. 12. The 
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stress data in these figures were derived from two types of strain gages: 
the ASTE gage, located above the shroud and close to the trailing edge; 
and the TIP gage, located close to the blade tip and about one third of 
a chord length beliind the leading edge. The strain gage located above 
the shroud at the maximum thickness point, and tlie gages below the 
shroud did noi register any appreciable stresses. 

The stress and displacement measurements are more directly re- 
lated to the blade motion than the pressure measurement. The results 
based on these two methcxls can, therefore, be used to evaluate the 
accuracy of results derived from pressure measurements. It can be 
observed from Figs. 11 and 12 that amplitude and phase distributions 
derived from pressure measurements are qualitatively correct despite 
a lack of precise quantitative agreement. However, even approximate 
knowledge regarding vibratory amplitude and phase distributions can be 
valuable. Such inlormation cannot be easily obtained from strain gage 
measurements because of their usually limited coverage. 

The agreement between stress and displacement amplitudes in 
Fig. 11 is very good except for blades 25 and 30. Anotlier unusual 
feature assoc laied with blade 25 is the significant difference between 
normalized ASTE and TIP stresses. Examination of Fig. 12 indicates 
that differences between stress phase-angles and displacement phase- 
angles are within 20'\ and that in particular, blades 25 and 30 are no 
exceptions. It is rioted that even if stress amplitude were in error, 
one would not exited the stress phase-angles to be in error. Examina- 
tion of Uie change in resistance of these strain gages subsequent to the 
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test run, however, indicates that according to the established practice 
the change was in&uificient to bring their accuracy into question |7|. 
Another cause for erroneous indication which cannot be ruled out is the 
suix?rsensitivity phenomenon |8]. Substantiation of tliis cause, however, 
would require further investigation. 

It appears, therefore, tliat {xisitive identification of the cause of the 
above inconsisiencies ca nnot be made at this iwint. However, the pro- 
ceeding disc ussion does iilustrate :hat in addition to providing displace- 
ments, the PES meihod can be usei in conjunction with stress measure- 
ments to provide some degree of redundancy and, therefore, increase 
the confidence in measurements. 

CONCLUSIONS 

1. Flutter can be detected i‘rom high-response aerodynamic data 
records. 

2. It is possible to obtain qualitatively correct distribution of 
vibratory amplitude and phase around the rotor from the data recorded 
by case-mounted static pressure transducers. 

3. Precise determination of flutter frequency can be made using 
two stationary transducers. 

4. Spectral displacement -amplitudes and phase angles during 
flutter can be obtained from the optical blade-tip ijosition measurements. 
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APPENDIX 

The S|X?ctral results presented in the main body of the paper were 
obtained bused on S()ecifications in Table 2. 

The factor of 16 in columns 3 and 4 is present because the tape 
siJeed on the play-back was reduced 16 times. The engine order fre- 
quency was determined indejxjndently for each data frame by numeri- 
cally detecting the blade passing frequency and then dividing it by the 
number of blades. Blade passing frequency was a predominant signal 
and it could bt' clearly detected despite tlie fact that it was beyond the 
filter cut-off frequency. The filter roll-off rate was 48 db per octave. 

Spectral analysis corresponding to Figs. 8 and 9 was performed 
on a large computer using an algorithm develoiJed by Singleton [9], 

A convenient feature of this program is that the number of points in a 
data sample does not have to be a power of two but it can be arbitrary. 
The program is included in tlie IMSL (International Mathematical and 
Statistical Libraries, Inc.) library. 

The remainder of the Apiiendix is concerned with the determination 
of phase angles relative to the start of a data sample. 

If the Fourier transform is applied to a digital data sample that con- 
tains a nomntegral number of cycles at a particular frequency, the re- 
sults will be not one but several peaks clustered around that frequency. 

It is usually sufficiently accurate to determine frequency based on the 
highest i^eaks. However,, the phase angle shifts by approximately 180° 
between peaks on each side of the true frequency. This can be seen 
by applying a sinusoidal input to the Fourier transform [5]. For a data 
sample of N jjoints the result is 
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1 ^ 1 _8ilLliL «'(N-1)/N (2) 

N N sin(c ll/N) 

where € = (u’q - a’),' (2 II /NT), u/q is the frequency of the sinusoidal input 
and T the sampling interval. When w * u)q, the magnitude of the 
normalized Fourier transform (I N)F(jw) is one. However, if the high- 
est tx?ak occurs at Wj < Wq, and there is also a significant con- 
tribution at w = u> 2 * a >2 > Wq, 

< sin c I II A. 

N sin(c jlI/N) 

where Aj and A 2 are amplitudes corresiX)nding to and 0^2 
respectively. Fraction c j, can therefore be estimated from this 
expression either numerically or graphically by plotting the left hand 
side. Once c j is determined it follows from equation (2) that the 
phase angle at u; = is equal to the phase angle at w = Wj minus 
€ jlI(N - 1)/N. 
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TABLE 1. - PHASE RELATIONSHIPS BETWEEN 


TWO TRANSDUCERS 47. 4° APART 


j Frequency, Cross •|X)wcr-si>ectrum 

E j |)hase, 

i deg 

Phase based | 
on Nj, 
deg , 

3.45 122.0 

4.45 173.6 

5.45 -144.2 

6.43 , - 93.7 

8.44 1 - 3.3 

; 11.46 150.0 

1 12.45 -158.3 

123.0 
170.5 
-142. 1 
- 94.7 
0.0 
142. 1 
-170.5 


TABLE 2. - PARAMETERS EMPLOYED IN SPECTRAL ANALYSES 




Sjx'ctral 

results 

Number of 
data ]X)ints 

Sampling 

frequency 

Filler cut-off 
frequency 

Number of 
data samples 

Rotational 

frequency 

Table 1 

2048 

819.2 • 16 

200- 16 

16 

127.31 

Figs. 1, 

3, and 4 

4096 

819.2 • 16 

200- 16 

16 

136.58 

Fig. 2 

4096 

819.2 • 16 

200- 16 

16 

116.78 

Fig. 5 

8192 

512 ' 16 

200- 16 

1 

127.31 i 

Fig. 8 

1026 

38E 


1 

127.31 

Fig. 9 

L .. .. 

IE 


1 

127.31 1 

1 


Unless noted otherwise, all frequencies are in Hz. 
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Fig. 1 Power spectrum of the wake-probe transducer signal 
during flutter. 





Fig. 2 Power spectrum of the wake-probe transducer signal; no 
flutter condition. 
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Fig. 3 Power spectrum of the midchord pressure transducer 
signal. 



Fig. 4 Power spectrum of the leading-edge pressure transducer 
signal. 
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Fig. 5 Linur spectrum of the pressure transducer signal. 



Fig. 6 Vibratory amplitude distribution from static pressure data. 
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H 9 . 7 Phase angle distribution from static pressure data. 



Fig. 8 Relative trailing edge displacements during flutter; 1 unit ■ 
0 . 71 mm. 
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Fiq. 12 Phase angle Jistribution at flutter frequenq^. 


N ASA‘L«wit'Ca«n*l 


